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ABSTRACT 

RX J1713.7— 3946 is the most remarkable very-high-energy 7-ray supernova 
remnant which emits synchrotron X-rays without thermal features. We made 
a comparative study of CO, Hi and X-rays in order to better understand the 
relationship between the X-rays, and molecular and atomic gas. The results 
indicate that the X-rays are enhanced around the CO and Hi clumps on a pc 
scale but decrease inside the clumps on a 0.1 pc scale. Magnet ohydro dynamic 
numerical simulations of the shock interaction with molecular and atomic gas 
indicate that the interaction between the shock waves and the clumps excite 
turbulence which amplifies the magnetic field around the clumps. We suggest 
that the amplified magnetic field around the CO and Hi clumps enhances the 
synchrotron X-rays and possibly the acceleration of cosmic-ray electrons. 

Subject headings: cosmic rays - X-rays: ISM - ISM: supernova remnants - ISM: 
individual objects (RX J1713. 7-3946) - ISM: clouds 
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Introduction 



RX J1713.7— 3946 is one of the most prominent supernova remnants (SNRs) emitting 
highest energy radiation covering the 7-rays and X-rays. The SNR is located relatively 
close to the Galactic center at (I, b) = (347?3, — 0?5), where contamination by the 
foreground/background features is heavy at any wavelength. The SNR was not known as 
a SNR in the radio continuum em ission and was discovered in the X-rays by the ROSAT 



( iPfeffermann fc Aschenbach 



19961 ) . The VHE (very high energy) 7- rays from the SNR w ere 



discovered and mapped by th e atmo spheric Cherenkov telescopes (jEnomoto et al. 



Aharonian et al. 



2004. 



2006a 



b 



2002 



20071 ). In particular, the H.E.S.S. observations resolved that 



the 7-rays distribution is shell-like with a diameter of 1° by 0.1-degree diameter PSF (point 
spread function). Since the 7-rays are detected at an energy range higher that 10 TeV, the 
cosmic ray (CR) protons producing the 7-rays may reach an energy range above several 100 
TeV, if the hadronic mechanism is working. The SNR is therefore an important candidate 
where the acceleration of the highest energy CR protons is best tested in the Galaxy. It is 
noteworthy that the X-rays of the SNR are purely non-thermal synchrotron emission; there 
are only two SNRs except for RX J1713.7— 3946, which show such no n-thermal X-rays, RX 



J0852.0-4622 (Vela Jr.) and 



Slane et al. 



2001 



Tian et al. 



HESS J1731-347, known to date (e.g., 



20101 ). 



Koyama et al. 



1997 



The distance of th e SNR was under controversy around 2000 (e.g. 



1997 



Wang et al. 



1997 



Slane et al. 



Koyama et al 



19991 ). It is now well established that the SNR is 



located at 1 kpc from us as first dete rmined by the rad ial velocity of the associated CO 



molecular gas at —7 km s 1 in Vlsr (IFukui et al. 



200 3J). This distance is confirmed by 



subse q uent observation s in X -rays and CO (jCassam-Chena'i et al. 



2004; 



20051)- 



Moriguchi et al. 



Moriguchi et al. 



(120051 ) showed further details of the 12 CO 



confirmed the interacting molecular gas identified by 



Fukui et al. 



' J =1-0 ) distribution and 
(120031 ) . Most recently, 
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Sano et al.l (2010, hearafter paper I) presented a comparison with the X-rays for the western 
rim of the SNR in 12 CO(J=2-l, 4-3) transitions, where most prominent mm/sub-mm CO 
clumps are distributed. These authors found that cloud core named "peak C" is a site of 
low-mass star formation and is surrounded by X-rays at pc scales, lending a further support 
for the distance of 1 kpc. Also, a comparison between the distribution of the Suzaku 
X-rays and the CO indicates that the dense cloud core shows an anti-correlation with the 
X-rays at sub-pc scales. The result suggests that the magnetic field is amplified around the 
dense cloud cores, leading to bright X-rays. It is noteworthy that the density distribution is 
generally highly inhomogeneous in the SNR as shown by these CO observations and such 



inhomogeneity remains to be incorporated ful 



particles with the interstellar medium (ISM) (IFukui 



y in treating the interaction of high-energy 



2008 



Zirakashvili fc Aharonian 



2010), 



rather than simply assum ing the uniform ISM distribution as in the previous works (e.g., 



Ellison et al 



2010. 



2012) 



Aiming at better understanding results of X- and 7-ray observations and their 
relationship with the neutral ISM, we are promoting a compre hensive comparison among 



the X-ra ys, the 7-rays and the ISM in 7-ray SNRs (e.g., Paper I; 



Fukui et al. 



paper II; 



Yoshiike et al. 



2012 . hereafter 



20131 ). Paper II presented an analysis of the interstellar protons 



both in molecular and atomic forms and showed a good spatial correspondence between the 
clumpy ISM protons and the VHE 7-rays. This provides a crucial support for the hadronic 
origin of the 7-rays instead of the leptonic origin. Now, in RX J1713.7— 3946, the next step 
is to establish the connection between the synchrotron X-rays and the ISM distribution. 
This work is aimed at better understanding the SNR shock-cloud interaction and obtaining 
a clue to the origin of the distribution of the synchrotron X-rays in the SNR. In the present 
paper we show a comparison of the spatial distribution among 12 CO(J=2-l), Hi and X-rays 
over the whole SNR in order to clarify the relationship between cold dense gas and the 
high-energy electrons. We will publish a separate paper that deals with a detailed spectral 
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analysis of the Suzaku X-ray results compared with the interstellar gas distribution (Sano 
et al. 2013b in prep.). 

The paper is organized as follows. Section 2 gives the description of the datasets of 
CO, Hi and X-rays and Section 3 gives the analysis. Section 4 gives discussion and Section 
5 the conclusions. 

2. Observations and data reductions 

2.1. CO 

The 12 CO(J=2-l) data at 230 GHz were taken with the NANTEN2 4 m telescope 
in the period f rom A ugust to November in 2008, and were published in Papers I, II and 



Maxted et al 



(120121 ). The front end was a 4 K cooled double sideband (DSB) receiver 
and a typical system temperature was ~250 K in the single sideband (SSB) including the 
atmosphere toward the zenith. The telescope had an angular resolution (FWHM) of 90" at 
230 GHz. We used an acoustic optical spectrometer has 2048 channels with a bandwidth 
of 390 km s _1 and resolution per channel of 0.38 km s _1 . Observations were carried out 
in the on-the-fiy mode with an integration time of 1.0 s or 2.0 s per grid, and provided a 
Nyquist-sampled 30" grid dataset. The ambient temperature load was employed for the 



intensity calibration. The absolute intensity scale was estimated 



KL object [5 h 35 m 14!52; -5°22'28"2 (J2000)] f lSchneider et al 



}y observing the Ori 



19981 ) and the main beam 



efficiency, ?? m b, was estimated to be 0.83. The rms noise fluctuations with 1.0 s and 2.0 s 
integrations were better than 0.66 K and 0.51 K per channel, respectively. The pointing 
accuracy was estimated to be better than ~15" by observing Jupiter every two hours. The 
image was smoothed with a Gaussian kernel with FWHM of 60". The observed area is 
shown in the Figured! 
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2.2. Hi 



Th e 21 cm Hi spectral data wer e taken from the Southern Galactic Plane Survey 



(SGPS; 



McClure-Griffiths et all 



1 120051 ) with the Australia Telescope Compact Array (ATCA) 
combined with the 64 m Parkes Radio Telescope. The combined beam size and the grid 
spacing of the dataset are 2.2 arcmin and 40", respectively. Also, the velocity resolution 
and typical rms noise fluctuations were 0.82 km s _1 and 1.9 K, respectively. We applied the 
correction for the Hi self-absorption by following the previous analysis of the cool Hi gas 
without CO emission (see Section 3.3 of Paper II). 



2.3. X-rays 



We used Suzaku archive data of RX J1713.7— 3946 taken from Data Archives and 
Transmission System (DARTS at ISAS/JAXA). The observations performed 15 pointings 
to ward the main fe a tures and 2 OFF pointings of R X J1713.7— 3946 and were published 



by 



Takahashi et al 



( 120081 ) and 



Tanaka et al. 



( 120081 ) expect for the 4 pointings observed 



in 2010 February. Previous and current observations are summarized in Table dj and the 
FoV of each observation is shown in Figure [TJ Activ e detector systems ab oard the Suzaku 



satellite are the X-ray Imaging Spectrometer (XIS; 



Koyama et al. 



20071 ) and the Hard 



X-ray Detector (H XD). The XIS consists o f four CCD cameras placed at the foci of X-ray 



Telescopes (XRTs; 



Serlemitsos et al 



20071 ). We analyzed on 



paper. The spaced-row charge injection (SCI; iNakajima et al 



XIS data in the pres ent 



2008. 



Uchiyama et al 



2009) 



was used in the latter 4 pointings (see also Table [TJ . Unfortunately, XIS 2 was on closed 
access since 2006 November 9, possibly owing to a micrometeorite damage. XIS showed 
an anomaly in Segment A on 2009 June 23. Thus, for the data in the latter 4 pointings 
we used XIS (expect for Segment A), XIS 1 and XIS 3. We used "cleaned event files" 
processed and screened by versions 2.0 or 2.4 Suzaku pipeline depending on observation 
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dates. We performed data reduction with the version 6.11 of the HEAsoft tools. 

Figure [2] shows mosaic images of RX J1713.7— 3946 which were constructed by using 
the data from XIS 0+1+2+3. The left and right panels show soft band (1-5 keV) and hard 
band (5-10 keV) images, respectively. The unit for the images is 10 -4 counts s _1 pixel -1 , 
and the pixel size is ~16.7". We smoothed the images using a Gaussian kernel with a 
FWHM of 45". We subtracted the non X-ray background (NXB) using xisnxbgen, which 
estimates NXB count rate based on night Earth observation data. Then, we correctet for 



XRT vignetting effects by simulating flat field images with xissim (jlshisaki et al. 



20071). In 



the soft band image, locations of the two bright point-like sources, 1WGA J1714.4— 39 45 



and 



WGA J1713.4-3949, are indicated by white circles. 



Pfeffermann fc Aschenbach 



(119961 ) claimed association of the former source with a Wolf-Rayet star. The latter is 



thoug ht to be a neutron star because of its X-ray spectral characteristics (lLazendic et al 
20031 ). 



3. Analysis 



3.1. Large scale CO, Hi and X-ray distributions 

Figure [3] shows four overlays of the 12 CO(J=2-l) distribution and X-ray images in the 
energy bands 1-5 keV (a, b) and 5-10 keV (c, d), respectively. A Vlsr range of CO from 
—20.2 to —9.1 km s _1 is shown in Figure [3] (a) and (c) from —9.1 to 1.8 km s _1 in Figure 
[3] (b) and (d), respectively. These overlays show a trend that CO and X-rays h ave a good 



correl ation at a pc scale. Th e CO p eaks (hereafter CO clumps) are identified by 



(2003:) and 



Moriguchi et al. 



Fukui et al. 



(120051 ) except for the clumps Dw and Osw that are newly 
identified in the present work. It is remarkable that most of the X-ray features are found 
toward CO clumps with only two exceptions, one of which is toward the candidate stellar 



remnant at (a J20 oo, feooo) = (17 h 13 m 28 s , -39° 49' 48") and the other is a Wolf-Rayet 
star (aj2ooo; feooo) = (17 h 14 m 30 s , —39° 46' 00"). In addition, ~10 X-ray point souces are 
distributed in the pointings but no excess can be seen toward their directions in Figure [21 

We see that the correlation between the X-rays and CO is good at a scale of pc, but 
at finer scales of 0.1 pc the peak of a CO clump tends to be shifted from the X-ray peak. 
Except for the clumps C, L and Osw, which are located inside of the SNR boundary, 
most of the CO clumps (A, B, D, Dw, E, G, O, Q and R) are embedded on the outer 
boundary of the SNR shell. The observed properties of the CO clumps are listed in Table 
[2J column (1): clump name; (2)-(3): position of the observed point with the maximum of 
12 CO(J=2-l) intensity; (4)— (6): observed properties of the 12 CO(J=2-l) spectra (radiation 
temperature, Vr,sR and FWHM line wid th); (7)-(9): size, luminosity and mass of the clump 



(jMoriguchi et al. 



2005 



Sano et al. 



20101 ); (10) : associated IRAS point source. 



Figured] (a) shows the same XIS mosaic image (1-5 keV) as Figure [2] (a), but the color 
scale is changed to stress the region of the low photon counts. We find the diffuse X-ray 
emission from the whole SNR has a good signal to noise ratio thanks to the low background 
of Suzaku. Figured] (b) shows the enlarged view in the SE-rim overlaid with the Hi proton 
column density contours (white). The integration range is —20.0 — 11.0 km s -1 . The lowest 
contour level and the contour interval in the Hi proton column density are 2.0 x 10 21 cm -2 
and 0.1 x 10 21 cm -2 , respectively. We apply the Hi self-absorption by following a previous 
analysis in Paper II. It is remarkable that the distribution of the self-absorption corrected 
Hi is complimentary to the X-ray peaks in the low-photon-count region in the SE-rim. We 
define a straight line PQ which passes through the center of the SNR and the X-ray peak 
in the SE-rim (see Section |3~2]) . 
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3.2. Detailed Comparisons with the X-Rays 

In order to illustrate the observational trend in Section 13.11 we made a comparison 
between CO clumps and X-rays for twelve CO clumps and present the results in Figures [5j 
The crosses in each image indicate the center of gravity of the CO clumps listed in Table [3] 
columns (2)-(3), which is weighted by the CO integrated intensity, differs somewhat from 
the CO peak position in Table |2] columns (2)-(3). In the right panels, we plot radial profiles 
of the CO integrated intensity and X-ray counts for the 1-5 keV and 5-10 keV bands 
centered on the center of gravity. Superposed are the Suzaku X-ray images in 1-5 keV 
and 5-10 keV. The higher energy band is free of absorption and the consistency between 
the two energy bands indicates that the absorption is not significantly affecting the X-ray 
images. Each arrow indicates the direction of the center of the SNR. The dashed magenta 
circles represent radii of 0?06 and 0?12 around the center of gravity for each CO clump. 
Right panels in Figures [5] show the radial profiles around each molecular clump in the 
12 CO(J=2-l) integrated intensity and the Suzaku two energy bands. All these features, 
except for the clumps D and Dw, clearly show that the X-ray emission is enhanced on 
the surface of each CO clump. This picture is consistent with the previous result on the 
clump C (Paper I). The averaged behavior is that the CO clumps have a radius of 0.04 ± 
0.01 degrees and the X-rays are distributed with a separation 0.06 ± 0.03 degrees from the 
center of each clump (see also Table [3]). 

Also, we have made a comparison between Hi clumps and X-rays around the SE-rim 
and present the results in Figure |6j The crosses indicate the Hi column density peaks on a 
straight line PQ (Section l3.ip . Radial profiles are also plotted. The trend of the relationship 
between Hi and X-rays is similar to that between the CO and X-rays, since the Hi peak 
trends to be shifted from the X-ray peak. 

In Figure [7] top panels show the same XIS mosaic images (1-5 keV) as in Figured] (a) 
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and middle and lower panels show the azimuthal angular distribution of the X-rays around 
each CO and Hi clump. The lines connecting the CO or Hi clump positions (shown in 
Figures 0, El and Table [3]) and the center of the SNR (I, b) = (347?3, — 0?5) or (aj 2 oocb <5j2ooo) 
= (17 h 13 m 34 s , —39° 48' 17") are defined as the origins of the azimuth angles, which is 
measured counterclockwise. The azimuthal angular distribution of the X-rays around each 
CO/Hl clump is estimated with respect to the direction of the center of the SNR. In order 
to do this, we estimate the background level of the X-rays to be ~66 x 10~ 4 counts s _1 
arcmin -2 , an average of the minimum value of the 1-5 keV X-ray profiles in Figure [71 This 
value is consistent with the nearly flat X-ray distribution inside the SNR. By setting a level 
at 1.5 times this X-ray background tentatively, we measured the azimuthal angular extent 
of the X-rays in each clump by using the 1-5 keV band image, because the signal to noise 
ratio is better in the 1-5 keV than in the 5-10 keV band. The results are shown in Table [3] 
and Figure In most of the azimuth distributions of the X-rays around each clump X-rays 
are enhanced at the azimuthal angle degrees (clumps A, B, D, Dw, E, G, O, Q, R and 
SE-rim). The three clumps inside the SNR, clumps C, L and Osw, are fully surrounded by 
the X-rays, which show a peak toward the azimuthal angle —180 degrees (clumps C and 
Osw) or +120 degrees (clump L). These distributions indicate that the clumps inside of the 
SNR are surrounded by the enhanced X-rays, and that those on the border of the SNR have 
enhanced X-rays only in the inner part of the SNR shell. 

4. Discussion 

4.1. Molecular and atomic clumps interacting with RX J1713.7— 3946 

Paper I showed that the synchrotron X-ray emission is enhanced around the CO clumps 
in the northwest region of RX J1713.7— 3946. In the present paper, we have compared the 
X-rays and the molecular and atomic clumps over the whole SNR and presented a schematic 
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image of these results in Figure [9j This figure indicates that the X-rays are enhanced 
around "all the CO and dense Hi clumps" within or on the outer boundary of the SNR 
shell. We infer that the CO clumps C, L and Osw survived the SNR blast waves because of 
their high density, being now embedd ed within th e SNR . This picture is consistent with the 



MHD numerical simulations made by 



Inoue et al. 



(120121 ) that the SNR shock waves cannot 



penetrate into the dense clumps, and the interaction produced turbulent and enhanced 
magnetic field around the clumps. The other CO clumps, for instance clump A, are also 
X-ray brightened only toward its inner side (about ±90 degrees in Figure [8]), indicating the 
shock interaction only at the surface of those CO and Hi clumps toward the center of the 
SNR. We see similar distributions of the X-rays in the high energy band (5-10 keV) that is 
hardly absorbed. The absorption in the X-rays therefore does not explain the depression. 
The morphology in Figures [7] and [S] indicates unmistakably that the X-ray emission is 
enhanced on the surface of the CO and Hi clumps and is produced by the interaction 
between the SNR blast waves and each CO/Hl clump. These CO and Hi clumps, now 
interacting with the SNR blast waves, were formed over a timescale of Myr by the stellar 
winds of the OB star which experienced a SN explosion 1600 yr ago (e.g., Papers I and II). 
Thus, the density in the SNR interior is expected to be very low. The observational results 
also indicate that there is no dense gas left in the interior of the SNR (see also Figure [3]). 
According to the numerical simulations, which studied the interaction of the ISM with the 
strong stellar w ind from O-type sta r, the gas density inside the evacuated wind bubble is 



n oj 



0.01 cm 3 (jWeaver et al. 



19771 ). This v alue is significant 



from the upper limit of the thermal X-rays (ITakahashi et al. 



y les s than 2 cm 3 , obtained 



20081 ). Therefore, the present 



results are consistent with that thermal X-rays has not been detected toward the inner 
region of the SNR. 
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4.2. Mechanism for the X-ray enhancement 



The X-rays from the SNR RX J1713.7— 3946 is purely synchrotron, which is spatially 
correlated well with the CO and Hi clumps on a pc scale but anti-correlated on a 0.1-pc 
scale as described above. If the X-rays trace the shock-amplified magnetic field and 
CR electrons, this trend is explained naturally a s the results of the inter action between 



Inoue et al. 



(2009 



20121 ) carried out 



the molecular and atomic clouds and the SNR. 
MHD numerical simulations of the interaction between the shock waves and the highly 
inhomogeneous neutral gas created by thermal instability to model RX J1713.7— 3946. The 
shock waves in the SNR propagate to the ambient cavity with density of < 1 cm -3 and the 
CO and Hi clumps having density of ~ 10 2 -10 3 cm -3 . Then the shock velocity in the dense 

clump becomes V^ciump = V shiam bi e nt(n amb ient/™cium P ) - 5 , wnere tne ambient density n amb i C nt 
= 1 cm~ 3 and the clump densi ty n r ] umr = 1 2 -10 3 cm -3 , and the shock is stalled when 



it collides with a dense clump. 



Inoue et al. 



( 120121 ) calculated that the temperature of the 
shocked dense gas becomes much lower than the temperature in the post-shock diffuse gas 
and the thermal X-ray line emission is strongl y suppressed, after the shock waves passed 



in it (see also Section 4.3 of 



Inoue et al. 



20121 ) . The shock velocity difference between the 



dense CO/Hl clump and the inter-clump ambient medium will become a factor of ~10-30. 
This velocity difference induces turbulence, which leads to turbulent dynamo action. As a 
result, the magnetic field is amplified to as high as 1 mG, which is consistent with the fie ld 



strength derived from rapid time variation of the X-ray filaments ( lUchiyama et al 



20071 ) 



while an altern ative is that the flu ctuation in the field orientation may explain the rapid 



time variation ( IHelder et al 



20121 ) . The synchrotron flux integrated in the line of sight i s 



Rybicki fc Lightman 



propo rtional to B 15 if the spectral index of electrons p is 2.0 (e.g., 

19791 ). So, it is possible to enhance the X-ray radiation around the CO and Hi clumps. The 

observed power of the X-ray emission around the CO and Hi clumps is 5-10 times higher 



than the inner part of the SNR. Then, the magnetic field around the CO and Hi clumps is 



13 



estimated to be 3-5 times higher than in the interior of the SNR, if the difference is due 
to only the magnetic field amplification. Therefore, the averaged magnetic field around 
the C O and Hi clumps b ecomes 45-75 /zG if the initial strength is assumed to be 15 fiG 



(e.g. 



Tanaka et al. 



20081 ) . The average fi eld strength is also estimated by the widt h of 



synchrotron X-ray filaments as ~100 /zG ( iBell 



2004; 



Hiraga et al. 



2005 



Ballet 



2006J). Note 



that the dependence of the synchrotron flux on the magnetic field strength can be much 
more sensitive than the above-mentioned standard case, because the high-ene rgy electrons 



that contribute the X-ray synchrotron emission can be in the cut-off regime (IBykov et al 



20081 ). Moreover, such enhanced magnetic field 



the diffusive shock acceleration (DSA; e.g., 



Bell 



may l ead to more efficient accele ration in 



1978 



Blandford fc Qstriker 



19781 ) scheme. 



Finally, we discuss quantitive relationship between the CO or Hi clump mass and the 
X-ray enhancement. We estimate the clump mass (column (9) of Tabel[2]) interacting with 
the shock wave as defined by the total CO or Hi mass within the azimuth angle range of 
the X-rays (column (9) of Table [3]) for each clump and plot it in Figure [10] as a function of 
the X-ray peak intensity in Figure [71 Also, we estimate the mass of the SE-rim to be 81 
M on the assumption that it is a sphere of radius 1.5 pc. The result, shown in Figure [T0| 
indicates that the correlation between the clump mass and the X-ray intensity is good with 
a correlation coefficient of ~0.84, by excl uding the lower li mits for the CO mass of clumps 



D and O, for which mass is not given by 



Moriguchi et al 



( 120051 ). Therefore, the X-ray flux 



is considered to be roughly proportional to the interacting mass of each CO and Hi clump 
at pc scales. This result suggests that the ISM distribution is crucial in interpreting the 
non-thermal X-ray distribution. More case studies of distribution of the X-rays and the 
ISM will be desirable to confirm this property. 
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5. Conclusions 

We summarize the present work as follows; 

1. We have shown that all the major CO and Hi clumps with mass interacting with 
the shock wave in RX J1713.7— 3946 are associated with the non-thermal X-rays. 
The X-rays are enhanced within ~1 pc of the CO and Hi peaks, whereas at smaller 
scales down to 0.1 pc the CO peaks tend to be anti-correlated with the X-rays which 
decrease toward the CO and Hi clumps. We have shown a quantitive relationship 
between the CO or Hi clump mass interacting with the shock wave and the X-ray 
intensity. The result indicates a good correlation between the two. 



2. The present findings in 1) is compared with numerical simulations of M 



ID in a 



Inoue et al. 



( 12009 



2012|). 



realistic highly inhomogeneous density distribution by 
These simulations indicate that the magnetic field is amplified around dense molecular 
or atomic clumps as a result of enhanced turbulence induced by the shock interaction. 
We interpret that thus-amplified magnetic fields enhance the X-ray intensity, which 
depends on the 1.5-th power of the magnetic field strength. Such enhanced magnetic 
field may also lead to more efficient acceleration in the DSA scheme. 
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Table 1. Summary of the Suzaku archive data of RX J1713.7— 3946 



Pointing ID 



ObsID 



a J2000 


5j2000 


XIS Exp 


/h m sN 


(° ' ") 


(ks) 


17 12 17.0 


-39 56 11 


69 


17 11 51.5 


-39 31 13 


18 


17 12 38.0 


-39 40 14 


21 


17 12 38.2 


-39 22 15 


22 


17 11 04.5 


-39 40 10 


21 


17 11 05.1 


-39 22 10 


21 


17 14 11.6 


-39 40 14 


21 


17 14 11.4 


-39 22 15 


18 


17 14 11.8 


-39 58 14 


21 


17 12 17.6 


-39 18 50 


21 


17 15 44.5 


-39 40 10 


20 


17 11 50.8 


-39 31 00 


62 


17 13 14.0 


-40 14 22 


19 


17 12 39.8 


-40 01 50 


21 


17 15 39.0 


-40 00 47 


22 


17 09 31.9 


-38 49 24 


35 


17 09 05.1 


-41 02 07 


38 



Date 



SCI 



100026010 

1 501063010 

2 501064010 

3 501065010 

4 501066010 

5 501067010 

6 501068010 

7 501069010 

8 501070010 

9 501071010 

10 501072010 

11 504027010 

12 504028010 

13 504029010 

14 504030010 

OFF1 100026020 

OFF2 100026030 



2005 September 26 OFF 

2006 September 11 OFF 
2006 September 11 OFF 
2006 September 11 OFF 
2006 September 12 OFF 
2006 September 12 OFF 
2006 September 13 OFF 
2006 September 19 OFF 
2006 September 19 OFF 
2006 September 20 OFF 

2006 October 5 OFF 

2010 February 15 ON 

2010 February 16 ON 

2010 February 17 ON 

2010 February 17 ON 

2005 September 25 OFF 

2005 September 28 OFF 



Note. 



Th e deta ils o f pointing ID from to 10, OFF1 and OFF2 were also shown in 



Takahashi ct al 



(2008) and 



Tanaka et al 



(2008) 
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Table 2. Properties of 12 CO(J=2-l) Molecular Clumps 



Name 


"J2000 


<5.J2000 


T R * 


•'peak 


AVlsr 


Size 


Lco^ 


Mass 1 " 


IRAS'* 




(h m s"\ 


(° ' ") 


(K) 


(km s _1 ) 


(km s- 1 ) 


(pc) 


(10 2 K km s' 1 pc 2 ) 


(M@) 


point sources 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


A 


17 11 38.4 


-39 58 46.9 


6.6 


-10.0 


1.5 


1.6 


32.5 


686 


17082-3955 


B 


17 12 26.8 


-40 05 55.5 


3.3 


-8.1 


4.5 


1.1 


9.0 


190 




C 


17 12 27.0 


-39 54 58.0 


7.5 


-11.9 


4.6 


1.2 


18.8 


397 


17089-3951 


D 


17 11 32.5 


-39 30 03.9 


3.3 


-9.3 


1.8 


1.3 


13.9 


292 


17079-3926 


Dw 


17 11 34.3 


-39 32 31.2 


3.1 


-1.1 


6.0 


1.1 








E 


17 11 21.0 


-39 47 24.4 


4.3 


-12.3 


2.7 


1.6 


7.5 


159 




G 


17 10 56.3 


-39 45 20.6 


2.8 


-11.5 


4.8 


1.5 


14.6 


307 




L 


17 12 30.2 


-39 28 14.6 


3.2 


-11.7 


6.0 


1.1 


17.6 


370 




O 


17 13 46.0 


-39 26 28.6 


1.2/1.9 


-4.6/1.0 


4.7/3.8 


1.2 


2.9 


61 




Osw 


17 13 08.9 


-39 36 43.2 


3.5 


-1.3 


1.3 


1.1 








Q 


17 15 11.7 


-39 26 47.4 


2.8 


-2.2 


2.8 


1.1 


5.1 


108 




R 


17 15 32.0 


-39 39 28.5 


3.2 


-3.1 


2.2 


1.1 


3.2 


67 





Note. — Col. (1): Clump name. Cols. (2)-(3): Position of the observed point with the maximum of CO(J=2— 1) intensity. 
Cols. (4)— (6): Observed properties of the 12 CO(J=2— 1) spectra obtained at the peak positions of the CO clumps. Col. (4): Peak 
radiation temperature Tr» . Col. (5): Vpe^ derived from a single Gaussian fitting. Col. (6): FWHM line width AVlsr- Col. 
(7): Size defined as (A/it) ' 5 x 2, where A is the total clump surface area defined as the region surrounded by values of peak 
integrated intensity. If the contour is unclosed, the boundary is defined as the intensity minimum between the nearby peaks. Col. 
(8): The CO luminosity of the clump Leo estimated by 12 CO(J=1-0). Col. (9): Total mass of the clumps derived by using 
the relation between the molecular hydroge n column den sity 7V(H2) and the 12 CO(J=1-0) intensity W( 12 CO), JV(H2) = 2.0 X 



10 20 [H/( 12 CO)/(K km s - 1 )] (cm~ 2 ) jBertsch et al 



' Moriguchi et al 



(2005) 



Sano et al 



(2010) 



19931 ). Col. (10): Associated IRAS point source. 
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Table 3. Results of Radial and Azimuthal Distribution 



Name 


"J2000 


<5j2000 


Vlsr 


Radius 


Separation 


Peak Intensity 


Angle (Fraction) 


Clump Mass 




/h m s \ 


(° ' ") 


(km s" 


- 1 ) 


(degree) 


(degree) 




(degree). 


%) 


(M@) 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


(8) 




(9) 


A 


17 11 39.1 


-39 59 22.2 


-14.5- 


-5.7 


0.050 


0.065 


239.5±0.9 


-90- 90, 


50) 


313 


B 


17 12 25.3 


-40 05 37.5 


-12.7- 


-3.6 


0.050 


0.025 


140.9±0.6 


-120- 60, 


50) 


95 


C 


17 12 25.3 


-39 55 07.4 


-16.6- 


-7.4 


0.040 


0.035 


254.3±0.9 


-180-180, 


100) 


397 


D 


17 11 31.9 


-39 29 13.6 


-14.2- 


-4.6 


0.040 




288.1±0.9 


-120-150, 


75) 


219 


Dw 


17 11 12.8 


-39 32 43.5 


-7.2- 


1.8 


0.040 




274.2±0.9 


-120- 90, 


58) 




E 


17 11 22.2 


-39 47 34.0 


-15.1- 


-9.8 


0.050 


0.025 


170±2 


-150- 90, 


67) 


106 


G 


17 10 54.5 


-39 46 25.7 


-16.4- 


-6.8 


0.050 


0.075 


152.9±1.2 


-60- 30, 


25) 


77 


L 


17 12 28.2 


-39 28 45.9 


-17.7- 


-5.7 


0.045 


0.095 


257.8±0.9 


-180-180, 


100) 


370 


O 


17 13 49.7 


-39 26 54.9 


-9.8- 


1.5 


0.040 


0.095 


182±2 


-120- 90, 


58) 


36 


Osw.... 


17 13 22.4 


-39 35 50.6 


-2.7- 


0.1 


0.030 


0.085 


188±2 


-180-180, 


100) 




Q 


17 15 09.6 


-39 38 48.7 


-5.1- 


0.3 


0.040 


0.025 


140±2 


-120- 90, 


58) 


63 


R 


17 15 32.4 


-39 39 28.5 


-5.3— 


-1.2 


0.035 


0.025 


86±2 


-120- 90, 


58) 


39 


SE-rim 


17 16 09.3 


-40 03 10.6 


-20.0— 


11.0 






90.3±1.1 


-60- 60, 


33) 


27 



Note. — Col. (1): Clump name. Cols. (2)-(3): Position of the center of gravity with the 12 CO(J=2-l) integrated intensity (except 
for the SE-rim). Cols. (4): Integration range for the estimation of the 12 CO(J=2-l) or H I (only SE-rim). (5): Radius of molecular 
clump defined as the radial distance from the center of gravity of CO clumps to the point where the half values of peak integrated 
intensity. (6) Separation of the X-ray peak intensity from each CO clump. (7) X-ray peak intensity around each clump shown in Figure 
[3 (unit is 10 -4 cnt. s — 1 arcmin -2 ). (8) Azimuth angle extent of the X-rays (1—5 keV), which more than background level shown in 
the text. (9) Clump mass interacting with the Shockwave defined as the total CO or H I mass within the azimuth angle range of the 
X-rays (1-5 keV) for each clump. 



-39 1 



£ -40> 



-41 l 
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Fig. 1. — Suzaku FoV of each observations toward the RX J1713.7— 3 9 46 ov erlaid on the 



Suzaku XIS 0+2+3 mosaic image (1-5 keV) taken from iTanaka et al.l (120081 ). The small 



squares correspond to the FoV of the XIS. The numbers indicated in the XIS FoV are 
pointing IDs used throughout this paper (see also Tabled]). We also plotted observed area 
in the 12 CO(J=2-l) enclosed by the white dash-dotted lines. 
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XIS 0+1+2+3 (1-5 keV) 



-39°15 



3 -39° 45 



-40°15 



17 h 16 m 17 h 12 m 

Right Ascension (J2000) 



28.0 



21.0 




14.0 i2 o 




1 7 h 12 m 
Right Ascension (J2000) 



Fig. 2.— Suzaku XIS (XIS 0+1+2+3) mosaic images of RX J1713. 7-3946 in the energy 
bands (a) 1-5 keV and (b) 5-10 keV. The color scale indicates count rate on a linear scale. 
The color bar numbers are in units of 10~ 4 counts s" 1 pixel -1 with a pixel size of ~16.7". 
Both images are smoothed with a Gaussian of 45". The positions of the two point sources, 
1WGA J1714.4-3945 and 1WGA J1713.4-3949, are shown with white circles (see the text). 
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1/lsr: -20.2 — 9.1 kms 



-39°15 



3 -39°45 



-40°15 




1/lsr: -9.1-1.8 kms 



-39° 15' 



3 - 39° 45' 



17 h 16 m 17 h 12 m 

Right Ascension (J2000) 



-40° 15' 




21.0 



14.0 m 



17 h 16 m 17 h 12 m 

Right Ascension (J2000) 



1/lsr: -20.2 — 9.1 km s 



■39° 15' 




VLsr: -9.1-1 .8 kms 



■39° 15' 



- 39° 45' 



-40° 15' 



17 h 16 m l/huP 1 

Right Ascension (J2000) 




17 h 16 m 17^12" 

Right Ascension (J2000) 



Fig. 3. — Pair of 12 CO(J=2-l) velocity channel maps (white contours) superposed with 
the Suzaku XIS mosaic image in two energy bands (a, b: 1-5 keV and, c, d: 5-10 keV) in 
color scale. The velocity ranges are (a) (c) —20.2 — 9.1 km s^ 1 and (b) (d) —9.1-1.8 km 
s _1 , respectively. The lowest contour level and the contour interval of CO are 3.1 K km s _1 
(~ 3a) in (b) and (d). In (a) and (c), the contour level are 3.1, 6.2, 9.3, 12.4, 15.5, 21.7, 
27.9, 31.0 K km s _1 . The CO clumps, A-D, Dw, E, G, L, O, Osw, Q and R are discussed in 
Section |U 
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XIS 0+1+2+3 (1-5 keV) 



-39°15 



3 -39°45 



-40°15 




17 h 16 m 17 h 12 m 

Right Ascension (J2000) 



6.0 B 



1/lsr: -20.0—1 1 .0 km s 



•39° 30' 



I -39° 50' 



■40° 10' 




6.0 m 



17 h 17 m 



1 7 h 15 m 
Right Ascension (J2000) 



1 7 h 13 m 



Fig. 4. — (a) Same XIS mosaic image (1-5 keV) as Figure |2] (a), but the color scale is 
changed to stress the region of the low photon counts. The scale is linear and units are in 
10 -4 counts s -1 pixel -1 . The lowest contour level and the contour interval are 3 and 1 x 
10 -4 counts s -1 pixel -1 , respectively in gray scale. This image shows X-rays by a factor of 
~2 fainter than that in Figure [31 (b) Enlarged view in the SE-rim. The white contours 
indicate the distribution of Hi proton column density (self-absorption corrected; see also the 
text and Paper II). The lowest contour level and the contour interval in Hi proton column 
density are 2.0 and 0.1 x 10 21 cm -2 , respectively. The velocity renge is —20.0 — 11.0 km s _1 . 
The crosses indicate the Hi column density peaks on a straight line PQ passing through the 
center of the SNR and the X-ray peak in the SE-rim. 
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Clump A (Hsu: -14.5 - -5.7 km s' 1 ) 




(min.: 3.8 K km s"\ step: 5.6 K km s' 1 ) 



l7 h 12 m 
Right Ascension (J2000) 



17 h 11 m 




1.5 2.0 



-•- NANTEN2 l! CO(J=2-1) 
-B- Suzaku 1-5 keV 
-A- Suzaku 5-10 keV x 8 



-'_:-'. ' : -f h 



A A -A— A- 



ClumpB (I4sn:-12.7--3.6kms 



7 h 12 m I7 n 11 m 

Right Ascension (J2000) 

(min.: 2.8 K km s"\ step: 4.7 K km s" 1 ) 



0.00 0.02 0.04 0.06 0.0 
Radius [degree] 



0.10 0.12 



1.5 2.0 



-40°00' 



-40°10' 



-39°50' 



-40°00' 




17 h 13 m 



17 h 12 
Right Ascension (J2000) 



17 n 13 m 17 h 12 m 

Right Ascension (J2000) 



0.00 0.02 0.04 0.06 0.0 
Radius [degree] 



r lg. 5. Distribution of 12 CO(J=2-l) emission (white contours) superposed with the Suzaku 1-5 keV and 5-10 keV 

from the left. Each CO integration velocity range and contour levels are showed in the top of the left and the middle panels, 
respectively. Each arrow indicate the direction of the center of the SNR. The crosses show the position of the center of gravity 
for the each CO clump (see also table |3). The dashed magenta circles represent radii of 0?06 and 0?12 around the center of 
gravity for each CO clump. Right panels show the radial profiles around each moecular clump in the 12 CO( J=2— 1) integrated 
intensity and Suzaku two energy bands (1—5 keV and 5—10 keV, in units of 10 — 4 counts s — 1 arcmin - 2 ) in Figure[2] The radial 
profiles from the 5-10 keV band has been scaled such that it has the same area as the 1-5 keV profile (scaled by a factor 8), 
for direct comparisons. 
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Clump D (I/Uh: -14.2 --4.6 km s" 1 ) 



(min.: 3.9 K km s' 1 , step: 2.0 K km s' 1 ) 



-39°25' 




0.0 0.5 



1.5 2.0 



• NANTEN2"CO(J=2-1) 
■ Suzaku 1 -5 keV 
Suzaku5-10keVx8 



^-T^ 



"*- T>-B-_, 



-A- a n 



17 h 12 m I7 h 11 m 

Right Ascension (J2000) 

Clump Dw ( Vlsh: -7.2 -48km s ') 



17"12 m 17 h 11 m 

Right Ascension (J2000) 



0.00 0.02 0.04 0.06 0.08 10 
Radius [degree] 



3.18 14 




1.5 2.0 



P ■ 4 



Right Ascension (J2000) 
Clump E (Vlsh: -15.1 --9.8 km s' 



17 n 11 m 30 s 



I 0.0 
17 n 10 m 30 s 
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-B- 
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Radius (degree) 




Right Ascension (J2000) 



Right Ascension (J2000) 



0.00 0.02 0.04 0.06 0.08 0.10 
Radius (degree) 



Fig. 5. — continued. 
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Clump G (l/Lsn:-16.4--6.8kms' 1 ) 



-39°45' 



D -39°55' 




(min.: 3.9 K km s' 1 , step: 2.0 K km s"') 



17 h 11 m 30 s 
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Clump O (Hsn:-9.8-4.5 km s" 1 ) 




(min.: 3.6 K km s\ step: 2.4 K km s" 1 ) 
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Fig. 5. — continued. 
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ClumpOsw (Vlsr: -27-0.1 km s'') 
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Fig. 5. — continued. 
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Fig. 6. — Distribution of emission Hi proton column density (white contours) superposed 
with the Suzaku 1-5 keV and 5-10 keV from the left. Hi integration velocity range and 
contour levels are showed in the top of the left and the middle panels, respectively. Each 
arrow indicate the direction of the center of the SNR. The crosses show the position of 
(«J2ooo, £j2ooo) = (17 h 16 m 9.3 s , —40° 3' 10.6") (see also table[3]). The dashed magenta circles 
represent radii of 0?06 and 0?12 around the crosses. Right panels show the radial profiles 
around the crosses in Hi proton column density and Suzaku two energy bands (1-5 keV and 
5-10 keV, in units of 10 -4 counts s _1 arcmin~ 2 ) in Figure [2J The radial profiles from the 
5-10 keV band has been scaled such that it has the same area as the 1-5 keV profile (scaled 
by a factor 8), for direct comparisons. 
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Fig. 7. — top panels: (a) Same XIS mosaic image (1-5 keV) as Figure H] (a) overlaid the 
circles centerd toward each CO/Hl clump. The white circle (radius = 0?12) and white lines 
from the center of the SNR are used for obtain the azimuth distributions shown in the lower 
panel, middle and lower panels: Azimuthal distributions of Suzaku XIS 1-5 keV (red 
square) and 5-10 keV (green triangle; scaled by a factor 8) integrated in rings centerd each 
clump. The lines connecting the center of gravity of CO/Hl clumps and the center of the 
SNR (I, b) = (347.3, —0.5) are defined as the origins of the azimuth angles. Here positive 
rotation is counterclockwise (e.g., clump Q in top panel). The horizontal solid lines indicate 
a background level of the X-rays (see the text). 
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Fig. 8. — The plots of the degree of angular extent subtended by enhanced X-rays (1-5 keV 
band) for each CO/Hl clump (see also Table [3]). Each black dot indicates the angles that 
the X-ray peak intensities around the each molecular clump. 
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Fig. 9. — Schematic image of the distribution of the molecular (CO) clumps (open crosses), 
atomic (Hi) clump (circle) and the X-ray radiation (shaded partial and full circles) super- 
posed with the Suzaku 1-5 keV X-ray boundary line of the SNR (gray contours). The black 
open crosses (clumps C, L and Osw) indicate that is able to survive against the shock wave. 
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Fig. 10. — Correlation plot between the X-ray peak intensity derived by Figure H and the 
clump mass interacting with the shockwave defined as the CO or Hi clumps mass within the 
azimuth angle range of the X-rays (column (9) of Table EJ) for each clump. The plots are 
fitted by the method of least squares. As a result, correlation coefficient to be ~0.84. We 
excluded the clumps D and O in the fitting (see the text). 



